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Abstract 

An explicit construction is provided for embedding n positive eigenvalues in the spectrum 
of a Schrodinger operator on the half-line with a Dirichlet boundary condition at the origin. 
The resulting potential is of von Neumann-Wigner type, but can be real valued as well as 
complex valued. 


1 Introduction 


Since the seminal paper of von Neumann and Wigner m, Schrodinger operators with embedded 
positive eigenvalues have always played a special role in spectral and scattering theory. In 
particular, several attempts have been made for generalizing the original result and for studying 
the stability of the embedded eigenvalues. It is obviously impossible to mention all references 
dealing with these questions, but let us mention a few which are related to our investigations 
El El El EHH], as well as the books p] Chap. 4.4] and uni App. 2, Chap. XI.8]. Despite these 
numerous works it seems to the authors of the present note that there is still some room left for 
discussing the case of n distinct embedded eigenvalues, especially for complex valued potentials. 

Given n positive numbers > /i 2 > • • • > n n > 0, we propose a very simple construction of a 
von Neumann and Wigner’s type potential V such that the corresponding Schrodinger operator 


— -pj + V on M + with Dirichlet condition at the origin, admits the eigenvalues • 


J hw 


The potential V as well as the eigenfunctions Vj are explicitly constructed and do not rely on an 
implicit equation. In addition, a family of n parameters can still be chosen arbitrarily, asserting 
once more that these eigenvalues are stable under suitable modifications of the potential, see for 
example mm- We emphasize that depending on the choice of these parameters, the potential 
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V can be real and thus leads to a self-adjoint operator, or can be complex valued. Note however 
that the leading term of the potential is of the form 

4 n 

V(r) = —^2 Hj sin(2//j-r) + 0(r~ 2 ) as r —> oo, ( 1 ) 

r j— i 

which implies in particular that the dependence in these parameters take place only in the 
remainder term at infinity. In fact, in the main statement below we exhibit the second term of 
the expansion of the potential and its dependence on the mentioned parameters. 

The possibility of constructing a potential V of the form ([Tj) with n positive eigenvalues 
has been known for a long time, see for example mum- Despite this fact, it appears to 
the authors that an explicit solution has never been provided. Let us however emphasize that 
our inspiration came from the paper [9] and from the book [ 6 ], in the special case when all 
mentioned parameters are equal to 1 (or more precisely when the matrix A of Theorem [T| is 
equal to the identity matrix). Moses and Tuan based their example on the Gel’fand-Levitan 
theory for inverse problem, see [7]. By adapting their idea, one can construct directly a potential 
having several positive eigenvalues, but it is not necessary to use the theory of Gel’fand and 
Levitan. Let also mention that an alternative construction has been proposed in m which leads 
to the possibility of embedding a finite or an infinite number of positive eigenvalues, but even 
in the finite case, the resulting potential is not of the form of the one we exhibit. In addition, it 
seems that most if not all previous works were dealing with real valued potentials only. 

In the last part of this note, we mention that our main result, obtained for Schrodinger oper¬ 
ator on M + with a Dirichlet condition at the origin, also leads to a similar result for Schrodinger 
operators on M 3 with the spherically symmetric potential V(\ ■ |). We also show that such a 
simple construction can only take place in M 3 and not in any other dimension if d is different 
from 1 or 3. 

Let us finally mention that part of the results of the present note was announced in m- 


2 The main result 


For any p\ > /x 2 > ■ ■ ■ > p n > 0 and any r£l, let us set 

s(r) := \ sm(mr), sin(/i 2 r),..., sin(/i n r)) 
c(r) := \ cos(pir),cos(p 2 r), • • •, cos(/i n r)) 

pr 

9 ij(r):= / sm(mp) sm(njp)dp, 

Jo 

where *(•) means the transposed vector (column vector). Note that g\j = gji, and by taking the 
equality 

sin {mp) sin(jLtjp) = - ( cos ((m - Hj)p) - cos ((/ij + Pj)p)^) 
into account, one also infers that 


9ij ( r ) 


hij (r) 

^ + ha(r) 


2 


for i 7 ^ j 
for i = j 


(2) 


with 


for i ^ j 


' sin ((m - nj)r ) _ sin ((^ + p,j)r) 
2(/i* — Tj) 
sin( 2 /Xjr) 

4/ii 


Both expressions for will be useful later on. 

We also define the n x n hermitian matrix G := and set l n for the n x n identity 

matrix. In the sequel we write G°°([0, oo)) for smooth functions on (—e, oo), for some e > 0, but 
restricted to the subset [0, oo). Finally, we write C* for C \ {0}, and use the standard notation 
f for the derivative of a function / with respect to its variable. Our main result then reads : 


Theorem 1. Let A = Diag(ai, • • • , a n ) be a diagonal n x n matrix with aj G {z G C* | > 

0 }. Then, (A + G(r)) is invertible for any r > 0, and by setting 


v(r) : = — [A + G(r)) 1 s(r) and V(r’) := 2(^(r), 

3=1 


the following properties hold: 

(i) V € L°°(M + ) fl C°°([0, oo)) and satisfies for r —>• oo 

4 n g n 

= f*3 sin (. 2 H r ) + 3 ( S a J N sin(2 Hj T ) + W (r)) + 0(W 3 ), (3) 

j = 1 3=1 

with W the real valued function given by 

n 2 71 

bC(r) = sin 2 (MjO) + 2 XT sin(/^r) cos (mr). 

3=1 *,3=1 

For j G {l,...,n}, t/ie component Vj of v belongs to G°°([0,oo)) and satisfies \vj(r)\ < 
Const./or any r G R+. In addition, Vj is a solution of ( — ^ + V)vj = pf-Vj, where 
the Dirichlet realization of the operator — ^ + V is considered in L 2 (M + ). 

Remark 2. Let us emphasize that the essential spectrum [0, oo) of the Dirichlet realization 
of --^pz + y, os well as its embedded eigenvalues pfj, are independent of any choice of the 
parameters aj, as long as the invertibility condition of A+ G{r) is ensured. This property is 
quite remarkable since the Vj’s and hence the potential V depend on A, as shown in Note that 
in this expression, we have emphasized the linear dependence on the parameters aj in the second 
order term, but the remainder term also depend on them. Note also that V is real-valued, and 
thus — + V is self-adjoint, if aj > 0 for all j G {1,... , n}, but is complex-valued if 9?(aj) / 0 

for some j. This persistence of the embedded eigenvalues under lower order modifications of the 
potential is consistent with the results contained in mw- 

Remark 3. For information and as we shall see in ©, © and in m, each Vj satisfies the 
expansions for r —> 0 

y j (r) = -a/Vi r + 0(r 4 ), 
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and for r —> oo 


Vj(r ) = sin (pjr) + — {aj sin (njr) + Y h jt (r) sin (mr)} + 0(r 3 ). 

l=i 

In a similar way, the asymptotic expansions of vl (r) as r —>• oo can be obtained from (USD, and 
the asymptotic expansion ofv"(r ) as r —> oo follows readily from the relation v" = (— pp- + V)vj. 

The following proof is divided into several small pieces. We use the notation M n ( C) for 
the set of all n x n complex matrices, and the notation (£,£) := ]C " =1 fjQ for the usual scalar 
product in C n . 


Proof, a) Let us first show the invertibility of A + G(r) for any r > 0. For r = 0, it follows from 
the definition that (7(0) = 0, and then A is invertible because each a 3 A 0. Assume now that 
there exists £ £ C n \ {0} which belongs to the kernel of A + G(r) for some fixed r > 0. Then 
one has 

0 = <£, (A + G(r))f) = i(£, $S(A)Q + <£, [»(A) + G(r)]£>. (4) 

By assumption on aj, it follows that 5R(A) > 0, while for the second real term one has 


<e,G(r)£> 




n 

Y & sin (mp) 

i= 1 


2 

dp > 0 . 


One infers from these estimates that the real term in Q can not be equal to 0, leading thus to 
a contradiction. 

Since A + G(r) is invertible for any r > 0, one sets 

v(r) := — (A + G(r)) 1 s(r) (5) 

and readily infers that v £ C' oo ([0, oo); C n ). 

b) Let us observe that for any fixed r > 0 and by the mean value theorem, there exists 
6 = 6{r) £ ( 0 , 1 ) such that gij(r) = rs\n(pn0r) sin (pjOr), from which one deduces that 




sin (piOr) 


l^iOr 


sin (pjOr) 


Pj9r 


p,iPj9^ ’ PiPj. 


Consequently, from this estimate and from the explicit expression provided in Q the following 
properties hold: 


J|l n + 0 ( 1 ) as r — > oo _ J|(l n + 0 (r 3 )) as r — » oo 
\ A + 0(r 3 ) as r —>■ 0 \ A + 0(r 3 ) as r —>-0, 


and for the inverse of this matrix, one deduces that 


(ATG(r ))- 1 


f 2 l n + 0 (r 2 ) as r —> oo 
[ d _1 + 0 (r 3 ) asr—» 0 . 


( 6 ) 


As a consequence of these estimates, one infers from the definition of v(r) and from the 
estimate | sin(//,;r)| < mr that for any j £ { 1 ,..., n} and any r > 0 


T 

vAr) I < Const.-~ 

3 w 1 _ 1 + r 2 


(7) 
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It thus follows that Vj E L 2 (M + ). 

c) Let us now consider v', the derivative of v with respect to its variable. Since 

v\r) = [A + G(r)) 1 G'^)(A +G(r)) 1 s(r) - (A + G(r))~ 1 s'(r), ( 8 ) 

one deduces from ([ 6 ]) and from the definition of G(r) that |u' (r)| < Const. (l+r) -1 , for any r > 0 
and any j E {1,..., n}. Note that a simple consequence of this estimate is that V E L°°(M+). In 
addition, the regularity property of V can easily be deduced from the corresponding properties 
of the functions v and s. 

d) Let us now check that the equality ( — 4~? + V^v = M 2 v holds, with M the diagonal 
n x n matrix Diag(pi, p 2 ; ■ ■ ■, p n )• In order to compute the expression —v" + Vv, observe that 
from the initial relation (A + G(r))u(r) = —s(r) one infers that 

G'(r)v(r) + [A + G(r))v'(r) = —Mc(r) 

and that 

G"(r)v(r) + 2 G'(r)v'{r) + (A + G(r))v"(r) = M 2 s(r), 

or equivalently that 

(A + G(r))v"(r) = M 2 s{r) — 2 G'{r)v'(r) — G"(r)v{r). 

From these relations, one then deduces that 

[A + G(r )) (v''(r) — \Vv\(r) + M 2 v{r) S j 
=M 2 s(r ) — 2 G'(r)v'(r) — G"(r)v(r) — V (r) [A + G(r)) v(r) 

+ M 2 (A + G(r))v(r) + [G(r), M 2 ]v(r) 

= — 2 G'{r)v'(r) — G"(r)v(r) + V(r)s(r) + [G(r), M 2 ]v{r ), (9) 


where [•, •] is used for the commutator of two matrices. 

On the other hand, observe that G'(r ) = (sin(//jr) sin(/r J r ))"^._ 1 = s(?’) t s(r), where we have 
used the identification of C n with the matrices M n i(C). As a consequence, one infers that 

G'(r)v'{r ) = (s(r),v'(r)) s(r), (10) 


and that 

G"(r)v(r) = (s(r),v(r)) s'(r) + (s'(r), v(r)) s(r). (11) 

Let us also observe that ([G(r),M 2 ]) = (p 2 — g 2 )gij(r ). In addition, by taking the equalities 

p 2 sin(p^p) = — [sin(/Ltj•) // ](p) into account, one also observes that 


pr pr 

(jp - Pi )9ij(r) = / sin (p,p) p 2 sin (ftp) d p- sin(p*p) sin(pjp) dp 

Jo Jo 

1 r f r 

3in(Pip)( - [sin^O'Kp)) + / [sin^-)'](/=») [ sin (M.j •)'](/>) dp 

J 0 Jo 
ir r r 

( - [sin(p r )'](p)) sin(pjp) -/ [sin(^-)'](p) [sin(/z J -) , ](/ 0 )(lp 

J 0 JO 


= [ - sin(jLtj-) sin(pj-) / + sin(p r )' sin(pj-)] (r). 
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Consequently, one has obtained that 

[' G(r),M 2 ] = — s(r)V(r) + s , (r)*s(r). (12) 

By inserting now the equalities (fTUj) . (fill) and (fT 2 l) into © one infers that 

(^4 + G(r)) (v"(r) — [Vv\(r) + M 2 n(r)j 

= — 2(s(r),v'(r)} s(r) — (s(r), u(r)) s^r) — (s^r), n(r)) s(r) + V(r)s(r) 

— (s\r), n(r)) s(r) + (s(r), n(r)) s r (r) 

= — 2f(s(r),v'(r)) + (s / (?’), v(r))^ s(r) + V(r)s(r) 

= ~ 2 [<s(-)^(-)> / ] (0»W + V{r)s(r) 

=0 

since — 2[(s(-),n(-)) , ](r) = — 2^sin(^j-) Vj(-)^ (r) = — V(r). Finally, since A + G(r) is 

invertible, one infers from the previous computation that v"(r) — [Vv](r) + M 2 v(r ) = 0, or 
equivalently that —v"(r) + \Vv\(r) = M 2 v(r), as expected. 

e) It has been shown in the point b) that Vj £ F 2 (M + ) for any j £ {1 and the 

equality Vj( 0) = 0 clearly holds. In addition, it follows from the estimate obtained in c) that 
Vj £ L 2 (M + ). Finally, since V £ L°°(M + ), as pointed out in c), one infers from the equality 
— v" = M 2 v — Vv that v" £ L 2 (M + ) as well. Thus, Vj belongs to the Dirichlet realization of the 
Laplace operator on M + , and this concludes the proof of the second statement of the proposition. 

f) It only remains to derive the asymptotic expansion (J3J). For that purpose, let us set 
H := {hij)] l j_ i ■ and let us use the notation || • || for the norms on C n and on M n ( C). By taking 
Q into account one gets 

A + G(r) = ^(l n +*(A + H(r))). 

Since ||IL(r)|| < C for all r > 0 with a constant C independent of r, one deduces that there 
exists ro > 0 such that for any r > r q 

V + ff(r)||<l/2 . 

r 

From the Neumann series, one then infers that 

(^4 + G(r )) -1 = -l n — ^(^4 + H(r)) + 0(r~ 3 ) as r ->• oo, (13) 

and by taking (fT3l) and (HI) into account, it follows that 

v(r) = “ s ( r ) + + H{r)}s{r) + 0(r -3 ), (14) 

and that 

2 4 

v'(r) =- Mc(r) 4 - -j {s(r) (s(r) s(r) + AMc(r) + H(r)Mc(r )} + 0(r~ 3 ). (15) 
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Finally, by putting together these information one obtains 

V{r) = 2 (s(r),i/(r)) + 2 (s'(r),v(r)) 

= 2 ^s(r), i/(r)) + 2 (Mc(r),r(r)) 

4 8 

= —(s(r), Mc(r)) H—^-(s(r), (s(r) (s(r) s(r) + AMc(r ) + id(r)Mc(r)}) 

— -(Mc(r), s(r)) + -^(Mc(r), {^4 + Ff(r)}s(r)) + 0(r~ 3 ) 

8 8 

= — (s(r), Mc(r)) 4 -^-{(s(r), AMc(r)) + (Mc(r), 7ls(r)) } 

+ -^{ || S (0H 4 + ( s ( r )i H(r)Mc(r )) + (Mc(r), i7(r)s(r))} + 0(r^ 3 ) 

4 , n A 8 . n A 8 

= — ^ /tj sin ( 2 AL r ) + — X] sin ( 2 /L' r ) + —^( r ) + 0(r -3 ), 
r i=i r i=i 

which gives the expansion ([ 2 ]). □ 

By using the standard relation between the Dirichlet Laplacian on R + and the restriction 
of the Laplace operator —A on R 3 to spherically symmetric functions (see for example [2] 
Sec. 11.3]), the previous result easily leads to a similar statement on R 3 . 

Corollary 4. The operator —A + V(| • |), with domain the Sobolev space % 2 (M 3 ), admits n 
eigenfunctions uj satisfying (—A + V(| • |))uj = /ijuj. These eigenfunction Uj are given by 
Uj(x ) \= Vj{\x\)f\x\ for any x £ R 3 and with Vj defined in Theorem\T\ 

Let us finally show that this construction is valid in R 3 only. Indeed, if we consider the d- 
dimensional Laplacian acting on spherically symmetric functions of the form Uj{x ) = a{r)vj{r ) 
with a,Vj € ^“((O, oo)), x £ R d \ { 0 } and r = \x\, then we find that 

[A Uj]{x) = a" (r) + ^—^u'j(r) 

= a(r)vj(r) + [2a' (r) + —— -a(r)Jvj(r) + [a"(r) + ——- a'(r)}vj(r ). 

In order to use a relation of the form — v” + Vvj = y'jV-j- it is thus necessary to impose that 

2a! (r) + -— -a(r) = 0 and a"(r) + -— -a'(r) = 0 . 
r r 

The first equation has the unique solution (up to constants) given by a(r) = r( -d+1 )/ 2 , and by 
substituting this solution into the second equation one obtains 

/// ^ i 1 // \ (d— l)(d —3) _(d+ 3")/2 n 

a (r) H-a (r =---r v ^ 11 = 0 

r 4 

which means that d can only be equal to 3 (the case d = 1 and a = const, clearly corresponds 
to even functions on R). 
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